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Abstract 

Kalman filters are often used to estimate the state variables of a dynamic system. 
However, in the application of Kalman filters some known signal information is often 
either ignored or dealt with heuristically. For instance, state variable constraints 
(which may be based on physical considerations) are often neglected because they do 
not fit easily into the structure of the Kalman filter. This paper develops an analytic 
method of incorporating state variable inequality constraints in the Kalman filter. 
The resultant filter truncates the PDF (probability density function) of the Kalman 
filter estimate at the known constraints and then computes the constrained filter 
estimate as the mean of the truncated PDF. The incorporation of state variable 
constraints increases the computational effort of the filter but significantly improves 
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its estimation accuracy. The improvement is demonstrated via simulation results 
obtained from a turbofan engine model. The turbofan engine model contains 3 state 
variables, 11 measurements, and 10 component health parameters. It is also shown 
that the truncated Kalman filter may be a more accurate way of incorporating 
inequality constraints than other constrained filters (e.g., the projection approach 
to constrained filtering). 

Key Words - Kalman Filter, State Constraints, Estimation, Probability Density 
Function, Gas Turbine Engines. 

1 Introduction 

For linear dynamic systems with white process and measurement noise, the Kalman 
filter is known to be an optimal estimator. However, in the application of Kalman 
filters there is often known model or signal information that is either ignored or dealt 
with heuristically [13]. This has resulted in recent efforts to incorporate constraints 
in the Kalman filter. For example, a projection method can be used to find the op- 
timal way to incorporate hard inequality constraints on the states [20, 21], Another 
way of incorporating constraints is to use a regularization method to enforce a soft 
limit on the changes of the state variables with respect to time [22], Yet another 
approach is the use of ridge regression to bias estimates with low certainty toward 
their constraints [5]. 

This paper presents a way to generalize the Kalman filter in such a way that 
known inequality constraints among the state variables are satisfied by the state 
variable estimates. The constraints that are imposed are hard constraints in that 
they are strictly enforced. However, in contrast to the projection method of con- 
straint enforcement [20, 21], the state estimates are not projected onto the constraint 
surface. Rather, the PDF that is computed by the Kalman filter is truncated at the 
constraint edges, and the constrained state estimate becomes equal to the mean of 
the truncated PDF. This idea is based on a previously published method [18] but 
has been modified to handle two-sided inequality constraints. 

The application considered in this paper is turbofan engine health parameter 
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estimation [6] . The performance of gas turbine engines deteriorates over time. This 
deterioration reduces the fuel economy of the engine. Airlines periodically collect 
engine data in order to evaluate the health of the engine and its components. The 
health evaluation is then used to determine maintenance schedules. Reliable health 
evaluations are used to anticipate future maintenance needs. This offers the benefits 
of improved safety and reduced operating costs. The money-saving potential of such 
health evaluations is substantial, but only if the evaluations are reliable. The data 
used to perform health evaluations are typically collected during flight and later 
transferred to ground-based computers for post-flight analysis. Data are collected 
each flight at the same engine operating point and corrected to account for variability 
in ambient conditions. Typically, data are collected for a period of about 3 seconds 
at a rate of about 10 Hz. Various algorithms have been proposed to estimate engine 
health parameters, such as weighted least squares [7], expert systems [4], Kalman 
filters [25], neural networks [25], and genetic algorithms [11]. 

This paper develops the truncation method of constrained Kalman filtering, and 
then applies it to the estimation of engine component efficiencies and flow capacities. 
Engine component efficiencies and flow capacities are referred to as health parame- 
ters. We can use our knowledge of the physics of the turbofan engine in order to 
obtain a dynamic model [2, 24], The health parameters that we try to estimate 
can be modeled as slowly varying biases. The state vector of the dynamic model 
is augmented to include the health parameters, which are then estimated with a 
Kalman filter [8]. We use heuristic knowledge of the health parameter dynamics 
to constrain their estimate. For example, we know that health parameters never 
improve. Engine health always degrades over time, and we can incorporate this in- 
formation into state constraints to improve our health parameter estimation. (This 
is assuming that no maintenance or engine overhaul is performed.) It should be em- 
phasized that in this paper we are confining the problem to the estimation of engine 
health parameters in the presence of degradation only. There are specific engine 
cases that can result in abrupt shifts in filter estimates, possibly even indicating an 
apparent improvement in some engine components. An actual engine performance 
monitoring system would need to include additional logic to detect and isolate such 
faults. 
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Section 2 derives the constrained Kalman filter. Section 3 discusses the problem 
of turbofan health parameter estimation, along with the dynamic model that we 
use in our simulation experiments. Although the health parameters are not state 
variables of the model, it is shown how the dynamic model can be augmented in 
such a way that a Kalman filter can estimate the health parameters [8, 12]. We 
then show how this problem can be expressed in a way that is compatible with 
the constraints discussed in the earlier section. Section 4 presents some simulation 
results based on a turbofan model linearized around a known operating point. We 
show that the truncated Kalman filter can estimate health parameters better than 
the unconstrained filter, and it can also estimate health parameters better than other 
constrained filters. Section 5 presents some concluding remarks and suggestions for 
further work. 

2 Constrained Kalman Filtering 

Consider the discrete linear time-invariant system given by 

x(k + 1) = Ax(k) + w(k) (1) 

y(k) = Cx(k ) + e(k) 

where k is the time index, x is the state vector, and y is the measurement. The sig- 
nals {u>(fc)} and {e(fc)} are uncorrelated zero mean Gaussian noise input sequences 
with covariances 

E[w{k)w T (rn)\ = Q6 km 
E[e(k)e T (m )] = R6 km 
E[w(k)e r (m)\ = 0 

where E[-] is the expectation operator and 8 km is the Kronecker delta function 
(8 km = 1 if k = m, 8 km = 0 otherwise). The Kalman filter equations are given as 
follows [1], 

K[k) = AZ(k)C T (CZ(k)C T + R)- 1 (2) 
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x(k + 1) = Ax(k) + Bu{k) + K{k){y{k) — Cx(k )) 
E(A + 1) = (AS(jfc) - K{k)C^(k))A T + Q 


where the filter is initialized with x(0) = i£[x(0)], and E(0) = E[{x — x(0))(x — 
x(0)) T ]. The Kalman filter estimate x(k) is a Gaussian random variable with a 
mean of x(k) and a covariance matrix of E(fc). 

Now suppose that we are given the s scalar constraints 

ai(k) < (f)J (k)x(k) < bi(k) i = 1, . . . , s (3) 

where Oj(fc) < bk(k). This is a two sided constraint on some linear function of the 
state. If we have only a one sided constraint, then we set a t (k) = — oo or bi(k ) = oo. 
Now suppose at time k that we have some estimate x(k) with covariance S(fc). The 
problem is to truncate the Gaussian PDF N(x(k),T,(k)) at the s constraints given 
in (3), and then find the mean x(k) and covariance S(fc) of the truncated PDF. 
These new quantities, x(k) and S(fc), become the constrained state estimate and its 
covariance. 

In order to make the problem tractable, we will define Xi{k) as the state estimate 
after the first i constraints of (3) have been enforced, and Ej(fc) as the covariance 
of Xi(k). We therefore initialize 

i = 0 (4) 

Xi{k) = x(k) 

%{k) = >:(/.•) 

Now perform the following transformation. 

Zi {k) = RW~ l / 2 T T {x(k) - Xi(k)) ( 5 ) 

where T and W are obtained from the Jordan canonical decomposition of Sj(fc). 

TWT t = ti(k) (6) 

We see that T is orthogonal and W is diagonal (therefore its square root is very easy 
to compute). Note that Zi{k) has a mean of 0 and covariance matrix of identity. 
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Next we use Gram-Schmidt orthogonalization to find the orthogonal R that satisifes 
RW^ 2 T T Mk) = [ 0 • • • 0 ] (7) 

With these definitions we see that the upper bound (3) is transformed as 


4>J {k)x{k) 

$ ( k)TW l / 2 R T Zi {k ) + $ ( k)xi(k ) 
Q i>J{k)TW l / 2 R T ) Zi {k ) 
{4>J (k)%(k)4>i(k )) 1 / 2 

1 0 • • • 0 1 Zi(k) 


bi[k ) (8) 

bi(k) 

bj{k) - (j)j{k) T Xj{k) 

(tfik^iWMk )) 1 / 2 

bilk ) - 4>i(k) T Xj(k) 

14? I fc)Si(fc)^(fc)) 1/2 

di(*0 


where dj(fe) is defined by the above equation. Similarly we can see that 

f 1 0 ... 0 1 z,(k) > a 4^ - 

> Cilk) 


where Cj(fc) is defined by the above equation. We therefore have the normalized 
scalar constraint 

Cilk) < 1 0 ••• 0 Zilk) < di{k) (9) 

Since z^k) has a covariance of identity, its elements are statistically independent of 
each other. Only the first element of ^(fc) is constrained, so the PDF truncation 
reduces to a one dimensional PDF trunction. The first element of Zj{k) is distributed 
as iV(0, 1) (before constraint enforcement), but the constraint says that Zilk) must 
he between Cj(fc) and dj(fc). We therefore remove that part of the Gaussian PDF 
that is outside of the constraints and compute the area of the remaining portion of 
the PDF as 

rdi(k) 1 1 r i 

/ — /== exp(— £ 2 /2) d( m - erf(dj(fc)/ V2) — erf(cj(A:)/ V2) (10) 

Ja(k) v 2tt 2 L J 

where erf(-) is the error function, defined as 

erf(/) -- [ exp(— t 2 /2)dt (11) 

V 71 " Jo 
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We normalize the truncated PDF so that is has an area of one, and we find that the 
truncated PDF (i.e. , the constrained PDF of the first element of Zi(k)) is given by 


PDF(C) = aexp(-C 2 /2) (12) 

V2 

a = f T 

0r eri(di(k)/V2) - erf(cj(fe)/\/2) 

We can compute the mean and variance of Zi(k ) as 

M = E[zi(k)\ (13) 

M (*0 „ 

= a Cexp(— C / 2 )d( 

Jci(k) 

= a exp(— cf (fc)/2) — exp(— d 2 (k)/2) 
cr 2 = E (zi(k) - n) 2 

rdi(k) 

= a (C — ^) 2 exp(— C 2 /2) d( 

Jci(k) 

= a exp(^c^(fc)/2)(c — 2/x) — exp(— dj(k)/2)(d — 2/j,) + /.i 2 + 1 

The mean and variance of the transformed state estimate, after enforcement of the 
first constraint, are therefore given as 

Zi+i{k) = n 0 ••• 0 (14) 

Cov(zi+i(k)) = diag(cr 2 , 1, • • • , 1) 

We then take the inverse of the transformation (5) to find the mean and variance of 
the state estimate after enforcement of the first constraint. 

x i+ i(k) = TW l / 2 R T z i+ i{k) + Xi{k) (15) 

t i+l {k) = TW l / 2 R T Cow{z i+l {k))RW l / 2 T T 

We then increment i by one and repeat the process of (5)— (15) to obtain the state 
estimate after enforcement of the next constraint. After going through this process 
s times (once for each constraint) we have the final constrained state estimate and 
covariance at time k. 

x(k ) = x s {k) (16) 

S(fc) = t s (k) 
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Figure 1 shows an example of a one-dimensional state estimate before and after 
truncation. Before truncation the state estimate is outside of the state constraints. 
After truncation, the state estimate is set equal to the mean of the truncated PDF. 
Figure 2 shows another example. In this case the unconstrained state estimate 
is inside the state constraints. However, truncation changes the PDF and so the 
constrained state estimate changes to the mean of the truncated PDF. 



— I I — 

-4 -2 0 2 4 


x 

Figure 1: The unconstrained estimate violates the constraints. The con- 
strained estimate is the centroid of the truncated PDF. 


3 Turbofan Engine Health Monitoring 

Figure 3 shows a schematic representation of a turbofan engine [16]. A single inlet 
supplies airflow to the fan. Air leaving the fan separates into two streams: one 
stream passes through the engine core, and the other stream passes through the 
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x 

Figure 2: The unconstrained estimate satisfies the constraints. Nevertheless, 
the truncation approach to constrained estimation shifts the estimate to the 
centroid of the truncated PDF. 

annular bypass duct. The fan is driven by the low pressure turbine. The air passing 
through the engine core moves through the compressor, which is driven by the high 
pressure turbine. Fuel is injected in the main combustor and burned to produce 
hot gas for driving the turbines. The two air streams combine in the augmentor 
duct, where additional fuel is added to further increase the air temperature. The air 
leaves the augmentor through the nozzle, which has a variable cross section area. 

Various turbofan simulation packages have been developed over the years [2, 3, 
10, 15]. The simulation used in this paper is a gas turbine engine simulation soft- 
ware package called MAPSS (Modular Aero Propulsion System Simulation) [16]. 
MAPSS is written using Matlab Simulink. The MAPSS engine model is based on 
a low frequency, transient, performance model of a high-pressure ratio, dual-spool, 
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Figure 3: Schematic representation of a turbofan engine. 


low-bypass, military-type, variable cycle, turbofan engine with a digital controller. 
The controller update rate is 50 Hz, and the component level model simulates the dy- 
namics of the engine components at a rate of 2500 Hz. The three state variables used 
in MAPSS are low-pressure rotor speed (XNL), high-pressure rotor speed (XNH), 
and the average hot section metal temperature (TMPC) (measured from aft of the 
combustor to the high pressure turbine). The discretized time invariant equations 
that model the turbofan engine can be summarized as follows. 

x(k + 1) = f[x(k),u(k),p(k)\ + w x (k) (17) 

p(k + 1) = p(k) + w p (k ) 

y{k) = g[x(k),u(k),p(k)\ + e(k) 

where k is the time index, x is the 3-element state vector, u is the 3-element con- 
trol vector, p is the 10-element health parameter vector, and y is the 11-elenrent 
measurement vector. The health parameters change slowly over time. Between 
measurement times their deviations can be approximated by the zero mean noise 
iUp(k). The noise term w x (k) represents inaccuracies in the system model, and e(k) 
represents measurement noise. An extended Kalman filter can be used with (17) to 
estimate the state vector x and the health parameter vector p. 
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The states, controls, health parameters, and measurements are summarized in 
Tables 1-4, along with their values at the nominal operating point considered in 
this paper (a power lever angle of 21° at zero speed at sea level). Table 4 also shows 
typical signal-to-noise ratios for the measurements, based on NASA experience and 
previously published data [14]. Sensor dynamics are assumed to be high enough 
bandwidth that they can be ignored in the dynamic equations. In Tables 1-4, LPT 
is used for Low Pressure Turbine , HPT is used for High Pressure Turbine , LPC is 
used for Low Pressure Compressor , and HPC is used for High Pressure Compressor. 


State 

Nominal Value 

LPT Rotor Speed 
HPT Rotor Speed 

Average Hot Section Metal Temperature 

7 264 RPM 
12 152 RPM 
1 533 °R 


Table 1: MAPSS turbofan model states and nominal values. 


Control 

Nominal Value 

Main Burner Fuel Flow 

Variable Nozzle Area 

Rear Bypass Door Variable Area 

2 454 lbm / hr 
343 in 2 
154 in 2 


Table 2: MAPSS turbofan model controls and nominal values. 

Constraints can be incorporated in the state estimator by using heuristic knowl- 
edge of the behavior of the health parameters. For example, it is known that health 
parameters never improve with time. It is also known that they degrade within a 
specific envelope. 

Pm(k) < p™ ax (T + 1), m € [1 - 10] (18) 

Pm(k ) > p™ m (fc + 1) 

This envelope constraint is in the linear form required in the constrained filtering 
problem statement (3) and is therefore amenable to the approach presented in this 
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Health Parameter 

Nominal Value 

Fan airflow 

1 

Fan efficiency 

1 

Booster tip airflow 

1 

Booster tip efficiency* 

1 

Booster hub airflow 

1 

Booster hub efficiency 

1 

High pressure turbine airflow 

1 

High pressure turbine efficiency 

1 

Low pressure turbine airflow 

1 

Low pressure turbine efficiency 

1 


Table 3: MAPSS turbofan model health parameters and nominal values. 
(*) The fourth health parameter is not yet implemented in MAPSS. 


paper. Note that this does not take into account the possibility of abrupt changes in 
health parameters due to discrete damage events. That possibility must be addressed 
by some other means (e.g., residual checking [6]) in conjuction with the methods 
presented in this paper. 

4 Simulation Results 

We simulated the methods discussed in this paper using Matlab. We measured a 
steady state 3 second burst of engine data at 10 Hz during each flight. Each of 
these routine data collections was performed at the single operating point shown in 
Tables 1-4, except the engine’s health parameters deteriorated a small amount each 
flight. The signal-to-noise ratios were determined on the basis of NASA experience 
and previously published data [14] and are shown in Table 4. The models on which 
this work was based are fairly comprehensive, so we assumed that the process noise 
for each component of the state derivative equation (17) was zero. However, in the 
Kalman filter we used a one-sigma state process noise equal to 0.005% of the nominal 
state values to allow the filter to be responsive to changes in the state variables. We 
also set the one sigma process noise for each component of the health parameter 


NASA/TM— 2006-2 14129 


12 




Measurement 

Nominal Value 

SNR 

LPT exit pressure 

19.33 psia 

100 

LPT exit temperature 

1394 °R 

100 

Percent low pressure spool rotor speed 

63.47% 

150 

HPC inlet temperature 

580.8 °R 

100 

HPC exit temperature 

965.1 °R 

200 

Bypass duct pressure 

20.66 psia 

100 

Fan exit pressure 

17.78 psia 

200 

Booster inlet pressure 

20.19 psia 

200 

HPC exit pressure 

85.06 psia 

100 

Core rotor speed 

12 152 RPM 

150 

LPT blade temperature 

1179 °R 

70 


Table 4: MAPSS turbofan model measurements, nominal values, and signal- 
to-noise ratios. 

to 0.01% of the nominal parameter value. These values were obtained by tuning. 
They were small enough to give reasonably smooth estimates, and large enough to 
allow the filter to track slowly time- varying parameters. In the enforcement of the 
constraints in (18) we chose the constraint envelope as follows. 

1. For the turbine airflow health parameters ( m € [7,9]), whose values increase 
with time (i.e., an increase corresponds to a degradation), p.™ ax (fc) was set 
equal to a linear-plus-exponential degradation that was initialized to zero (i.e., 
p™ ax (0) = 0) and reached a maximum of 6% after 500 flights, while p™ m (fc) 
was set equal to 0 for all k. 

2. For the other health parameters ( m € [1 — 6,8,10]), whose values decrease 
with time (i.e., a decrease corresponds to a degradation), p™ m (fc) was set 
equal to a linear-plus-exponential degradation that was initialized to zero (i.e., 
*C n (0) = 0) and reached a maximum magnitude of -6% after 500 flights, 
while p™ ax (fc) was set equal to 0 for all k. 

We simulated a linear-plus-exponential degradation of the 10 health parameters 
over 100 flights. The initial health parameter estimation errors were assumed to be 
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zero. The simulated health parameter degradations were representative of turbofan 
performance data reported in the literature [17]. 

Figure 4 shows a typical plot of the true deviation of health parameter 10, 
along with the constraint envelope, the unconstrained estimate, and the constrained 
estimate. It is seen that even though the unconstrained estimate lies within the 
constraint envelope, the constrained estimate is more accurate. Figure 5 shows a 
different type of example where the true health parameter deviation is closer to the 
constraint envelope. In this case there are times when the unconstrained estimate 
lies outside of the constraint envelope, but the enforcement of constraints forces the 
constrained estimate to remain within the envelope. 


Health Parameter #1 0 



Figure 4: In this example, constraint enforcement decreases the RMS estima- 
tion error from 12.2% to 9.2%. 

We ran 20 Monte Carlo simulations, each with a different noise history. We 
obtained estimates of the health parameters using three different methods. 
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Health Parameter #1 



Figure 5: In this example, constraint enforcement decreases the RMS estima- 
tion error from 13.4% to 6.6%. 

1. Unconstrained (standard) Kalman filtering. 

2. Constrained Kalman filtering using the projection approach [20, 21]. 

3. Constrained Kalman filtering using the projection approach and constraint 
tuning [23]. 

4. Constrained Kalman filtering using the truncation approach discussed in this 
paper. 

Table 5 shows the performance of the filters averaged over all 20 simulations. The 
standard Kalman filter estimates the health parameters to within 7.4% of their final 
degradations. The projection-based constrained filter estimates the health param- 
eters to within 6.5% of their final degradations. The projection-based constrained 
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filter with the addition of residual-based tuning estimates the health parameters to 
within 6.1% of their final degradations. Finally, the use of the truncation approach 
for constrained filtering estimates the parameters to within 6.1% of their final degra- 
dations. These numbers show the improvement that is possible with the truncation 
approach to constrained Kalman filtering. Although we may be able to get just as 
good performance using the tuned projection filter, a lot more tuning is required 
than with the truncation approach [23]. 



Estimation Error (%) 

Health Parameter 

Unconstrained 

Filter 

Projection 

Filter 

Tuned 

Filter 

Truncated 

Filter 

Fan airflow 

12.9 

9.2 

8.2 

7.5 

Booster hub airflow 

6.9 

6.2 

6.0 

5.8 

Booster tip airflow 

10.9 

10.6 

10.0 

10.5 

Booster tip efficiency* 

N/A 

N/A 

N/A 

N/A 

Booster hub airflow 

7.4 

6.8 

6.3 

6.4 

Booster hub efficiency 

3.8 

3.1 

3.0 

3.7 

High pressure turbine airflow 

4.3 

3.3 

3.2 

4.0 

High pressure turbine efficiency 

4.2 

3.8 

3.7 

4.2 

Low pressure turbine airflow 

3.6 

3.3 

3.2 

3.8 

Low pressure turbine efficiency 

11.3 

11.2 

11.1 

8.8 

Average 

7.4 

6.5 

6.1 

6.1 


Table 5: Health parameter estimation errors (percent) of the Kalman filters. 
(*) The fourth health parameter is not yet implemented in MAPSS. 


The improved performance of the constrained filter comes with a price, and that 
price is computational effort. The algorithm outlined in (5)— ( 15) requires Jordan 
decomposition and Gram-Schmidt orthogonalization. However, if the constraints 
of (3) are decoupled (as they are in our example) then the computational effort 
can be largely reduced by ignoring the cross-covariance terms in the state estimator 
and hence avoiding these matrix computations. In any case, computational effort is 
not a critical issue for turbofan health estimation since the filtering is performed on 
ground-based computers after each flight. 
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Note that the Kalman filter works well only if the assumed system model matches 
reality fairly closely. The method presented in this paper, by itself, will not work well 
if there are large sensor biases or hard faults due to severe component failures. A 
mission-critical implementation of a Kalman filter should always include some sort 
of additional residual check to verify the validity of the Kalman filter results [9], 
particularly for the application of turbofan engine health estimation considered in 
this paper [6]. 

5 Conclusion and Discussion 

We have presented a PDF truncation based method for incorporating constraints 
into a Kalman filter. If the system whose state variables are being estimated has 
known state variable constraints, then those constraints can be incorporated into 
the Kalman filter as shown in this paper. For the aircraft turbofan engine health 
estimation problem, the use of constraints generally improves the accuracy of health 
estimatation. At first this seems counterintuitive, since the unconstrained Kalman 
filter is by definition the minimum variance filter. However, we have changed the sys- 
tem by introducing state variable constraints. Therefore, the unconstrained Kalman 
filter is no longer the minimum variance filter, and we can do better with the con- 
strained Kalman filter. 

We have seen that the constrained filter requires more computational effort than 
the standard Kalman filter. This is due to the addition of s matrix decompositions 
that must be performed at each time step (one for each constraint). The engineer 
must therefore perform a tradeoff between computational effort and estimation ac- 
curacy. For real time applications the improved estimation accuracy may or may 
not be worth the increase in computational effort. 

The Kalman filter works well only if the assumed system model matches reality 
fairly closely. The constraint enforcement and constraint tuning methods presented 
in this paper will not work well if there are large sensor biases or hard faults due 
to severe component failures. A mission-critical implementation of a Kalman filter 
should always include some sort of residual check to verify the validity of the Kalman 
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filter results, particularly for the application of turbofan engine health estimation 
considered in this paper [6, 9]. 

Although we have considered only linear state constraints, it is not conceptually 
difficult to extend this paper to nonlinear constraints. If the state constraints are 
nonlinear they can be linearized as discussed in [19]. Further work could explore 
ways to optimally tune the constraints of the truncated Kalman filter. 
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